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1 .  Introduction.  In  this  paper  we  shall  investigate- sane 
general  properties  of  motions  of  a  compressible  fluid  under  the 
influence  of  a  gravitational  field.  Our  main  restrictive  hypotheses 
are  first  that  we  are  concerned  with -motions  (or  time  scales)  for 
which  the  dissipative  effects  of  viscosity  and  heat  condition  are 
not  important,  and  second  that  the  motions  are  small  deviations 
from  a  basic  state  of  rigidly  rotating  hydrostatic  equilibrium. 

We  also  assume- that  the  gravitational  field  is  externally  given, 
and  is  not  itself  affected  by  the  motion.  The  mathematical  model 
based  on  these  hypotheses  is  the  linear  theory  of  rotating  stratified 
non-dissipative  flow.  Perhaps  the  most  basic  mathematical  problem 
in  this  context  is  the  initial -boundary  value  problem:  the  values 
of  the  flow  variables  are  given  at  some  initial  time  in  a  region 
R,  and  their  subsequent  temporal  evolution  is  sought,  subject  to 
appropriate  boundary  conditions  on  the  boundary  B  of  R.  We  shall 
be  concerned  with  this  problem  in  the  case  of  the  boundary  condition 
corresponding  to  a  rigid  evelope  B,  i.e.  we  take  the  normal  component 
of  the  velocity  vector  to  be  zero  on  B.  As  usual  with  initial -boundary 
value  problems  for  linear  systems  of  partial  differential  equations 
with  time  independent  coefficients,  it  is  helpful  to  think  of  the 
problem  in  terms  of  the  "method  of  normal  modes , "  a  normal  mode  being 
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a  solution  of  the  differential  equations  which  satisfies  the  (homogeneous) 
boundary  conditions  and  varies  harmonically  with  time.  If  the  initial 
data  can  be  represented  as  a  superposition  of  the  initial  values  of 
the  various  normal  modes,  then  the  solution  to  the  initial -boundary 
value  problem  Is  simply  the  same  superposition  of  normal  mode  com¬ 
ponents,  each  evolving  in  time  at  its  own  frequency.  From  this 
point  of  view  there  are  two  parts  to  the  problem:  first  the  determin¬ 
ation  of  the  frequencies  and  the  spatial  dependence  of  the  various 
normal  modes,  and  second  the  representation  of  the  initial  data  in 
terms  of  their  initial  values.  In  many  problems,  for  instance  in 
heat  conduction,  the  first  part  is  by  far  the  more  difficult,  the 
representation  of  initial  data  being  readily  found  by  exploiting 
some  sort  of  spatial  orthogonality  of  normal  modes  of  different 
frequencies.  To  a  certain  extent  this  is  true  of  the  present  problem  - 
it  is  only  in  exceptional  cases  that  it  is  possible  to  describe  the 
oscillatory  normal  modes  explicitly  (indeed  even  their  existence 
is  not  an  easy  question),  though  they  do  have  a  kind  of  spatial 
orthogonality.  However  in  the  case  of  steady  flows,  normal  modes 
of  frequency  zero,  it  turns  out  to  be  possible  to  describe  them  quite 
generally.  But  in  the  present  case,  as  in  the  analogous  problem  for 
homogeneous  rotating  flow,  this  class  of  steady  flows  is  often  quite 
extensive  -  the  frequency  zero  is  highly  degenerate,  and  this  means 
that  the  second  part  of  the  problem  -  the  determination  of  exactly 
which  steady  component  is  excited  by  the  initial  data  -  is  by  no  means 
trivial.  Our  main  results  deal  with  this  aspect  of  the  problem.  After 
formulating  the  problem  and  describing  some  of  its  general  features, 
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we  shall  describe  wl 1  the  steady  flows.  Then  we  shall  establish  some 
basic  conservation  properties  possessed  by  all  flows  in  this  linear 
theory,  and  show  how  they  can  be  used  to  characterize  exactly,  from 
the  initial  data,  the  steady  component  of  the  solution  to  the  initial¬ 
boundary  value  problem.  We  conclude  with  some  examples  illustrating 
our  general  results  in  special  cases,  and  compare  the  present  problem 
with  its  analogue  for  rotating  homogeneous  flow. 

2.  Formulation .  The  basic  equations  of  our  mathematical  model 
are  those  of  a  compressible  fluid  without  •>  iscosity  or  heat  conduction, 
written  relative  to  a  system  rotating  about  the  z-axis  with  angular 
velocity  fi : 


u^  +  2fik  X  u  +  u  •  Vu  +  ^-Vp  +  Vcp  =  0  (1) 

Pt  +  V  •  (pu)  =  0  (2) 

s^.  +  u  •  Vs  =-•  0  (.5) 

P  =  f(p,s)  (4) 


Here  u,  P,  p  and  s  are  velocity,  density  pressure  and  specific 

2  2  2 

entropy  and  cp  =  -fi  (x  +  y  ;/2  +  is  the  "geopotential,"  cp^  being 
the  ordinary  gravitational  potential;  in  the  simplest  case  we  would 
have  cp^  =  gz,  but  we  do  not  require  any  specific  form  for  cp^,  nor 
need  we  specialize  the  equation  of  state  (4),  though  we  shall  assume 
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(as  is  appropriate  for  all  ordinary  fluids)  that  df/dp  >  0.  Let 

L  be  a  length  characterizing  the  size  of  the  regie®  R  containing 

the  fluid,  and  let  p  and  s  he  characteristic  density  and  entropy 

scales.  Then  introduce  dimensionless  variables,  indicated  by  asterisks, 

as  follows:  r  -  Lr^,  t  =  ft  t*,  u  =  HLu*.,  P  =  PP#,  s  =  ss# 

2  2~  2  2 

p  =  L  (1  pp#,  ?  =  !!  L  The  dimensionless  equations  are  then  the 

same  as  (1)  -  (4)  except  that  Q  is  replaced  by  1,  asterisks  are 
affixed  to  all  the  variables,  and  for  the  equation  of  state  we  write 
=  P  f(ss^,L  H  pp+).  Let  the  basic  state  about  which 
we  shall  linearize  have  u^.  =  0,  p^  =  p^,  p#  =  pQ,  s#  =  sQ;  the 
equations  then  show  that  pQ,  pQ  and  sQ  are  all  functions  of  9#, 
and  in  fact  PQ  =  -p’0(9#)  and  PQ  =  f#(sQ,p0) .  We  may  thus  consider 
the  basic  state  of  stratification  to  be  specified  by,  say,  P0(9„.), 
and  PQ  and  sQ  are  then  determined.  It  should  be  noted  that  while 
we  shall  usually  write  expressions  like  P  =  p  (9,. ) ,  this  should  not 

u  (J  * 

be  interpreted  quite  literally,  for  what  the  equations  lor  the  basic 
state  really  say  is  that  PQ  does  not  change  as  one  moves  along  a 
surface  of  constant  9*  in  the  flow  region.  But  the  constant  9K 

surface  may  be  cut  into  several  disconnected  components  by  the  flow 
region,  and  though  pQ  must  be  a  constant  within  each  component, 
the  different  components  may  have  different  values  of  this  constant. 

The  density  distributions  down  inside  each  of  two  potholes  in  the  bottom 
of  a  stratified  lake  need  not  agree.  The  linearized  equations  are 
now  obtained  by  setting  u*.  «  Eu. ,  p„  =  p_  +  ep., ,  p,,  =  p.  +  Ep, , 

“*X  U  X  ■*  0  X 

s„  =  s.  4  ts.  in  (1)  -  (4)  and  retaining  first  order  terms  in  E; 

*  (J  X 


e  may  be  regarded  as  a  sort  of  Rossby  number  characterizing  the 
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magnitude  of  the  initial  perturbation.  Suppressing  from  now  on  the 

subscripts  1  and  *,  and  writing  f  and  f  for  df„/&s  and 

s  p  ^  * 

df#/dp#  evaluated  at  the  basic  state,  these  linearized  equations  are: 

Uj.  +  2k  x  u  +  p'Vp  +  p^pVfp  =0  (5) 

Pt  +  V  -(PqU)  =  0  (6) 

st  +  u  •  VsQ  =  0  (7) 

P  =  fss  +  fpP-  (8) 

These  equations  are  to  be  considered  in  a  fixed  region  R,  with 
the  boundary  condition  that  u  •  n  =  0  on  the  boundary  B  of  R, 
n  being  the  unit  outward  normal  to  B.  For  the  initial  value 
problem  we  also  prescribe  u,  p  and  s  at  t  =  0. 

Hie  equations  (5)  -  (8)  possess  an  invariant  (i.e.  temporally 
conserved)  "energy  integral."  It  can  be  obtained  as  follows: 

&(ipoV)-  0^(0*-  Vt> 

=  p”1p(-P0V  •  U  -  p^u  ■  Vcp  +  s^fgu  •  vcp) 

=  -pV  •  u  -  p^pfpp^u  •  Vcp  =  -p  V  •  u  +  pfpU  •  Vcp, 

since  p!  =  f  s '  +  f  p!  from  the  original  equation  of  state.  But 
0  s  0  p  0 


we  also  have 
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±(  l  (gn-if  s2 

at  V  2  •  O'  a 


=  -(s')  ss!u  •  v<p  =  -sf  u  •  V<p, 
0  s  0—  s 


and  so 


a  (  i  -l.  2  i  ,  ,  -i. 

StU°OfPP  '2<S0'  f= 


s  )  =  -pV  •  u  +  Pu  •  V<P 


=  -V  •  (pu)  +  u  •  (Vfc  +  pW) 


-  (K)  -  f&fOoW*) 


using  (5).  Thus 


& |  {p0liil2  +  polfpp2  -  (s<y'lfss7  +  v  •  {p^  =  °- 


(9) 


Integrating  this  over  R  and  using  the  boundary  condition  we  find 
the  energy  conservation  formula: 


ac  s  JL  1  f 

at  at  2  JR 


7/ 

0  p 


(10) 


We  assume  that  f  >0  (density  increases  with  pressure  if  no 
P  “ 

heat  is  added),  and  we  shall  also  suppose  that  the  basic  state  is 

such  that  (s')  7  <  0,  i.e.  it  is  "statically  stable."  Usually 

u  s 

this  would  mean  s^  >  0  (specific  entropy  increases  with  altitude), 
since  most  ordinary  fluids  expand  on  being  heated  at  constant  pressure, 
i.e.  f  <  0.  With  these  assumptions  we  see  that  the  quadratic  form 
6  is  positive  definite  and  thus  the  perturbation  quantities  will 
remain  bounded  (in  norm  anyway)  if  they  are  so  initially.  It 
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is  of  course  to  achieve  this,  i.e.  to  rule  out  the  possibility  of 
convective  instability,  that  we  assume  static  stability  of  the  basic 
state. 

If  we  introduce  the  inner  product  corresponding  to  the  quadratic 
form  C  by 


<u1,p1,s1|u2,] 


*2,S2 


)=/EtPoUj 


<’o1yIP2  -  <si>'1V>2ldV  <u> 


(the  asterisk  denotes  the  complex  conjugate)  one  readily  verifies  that 
the  inner  product  of  two  (possibly  complex)  solutions  of  (5)  -  (8) 
is  also  invariant  in  time.  In  particular  this  implies  the  reality 
of  the  frequency  of  any  normal  mode  solution  of  (5)  -  (8),  and  the 
orthogonality,  in  the  sense  of  the  inner  product  (11),  of  any  two 
normal  modes  of  different  frequencies. 

We  note  also  that  the  energy  conservation  formula  (10)  implies 
the  uniqueness  of  the  solution  to  the  initial  value  problem,  for  the 
difference  of  two  soluti ons  with  the  same  initial  data  is  also  a 
solution  o'  the  equations  and  the  boundary  condition,  and  has  zero 
energy  initially.  Being  conserved,  its  energy  is  also  zero  there¬ 
after,  and  since  it  is  positive  definite  this  implies  that  the 
difference  flow  is  always  zero,  i.e.  the  original  solution  is 
unique.  A  proof  of  the  existence  of  a  solution  to  the  initial  value 
problem  appears  to  be  a  more  difficult  question  -  we  shall  merely 


assume  it  here. 
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2_.  Steady  flows.  The  time -independent  solutions  of  equations 
(5)  -  (8)  can  be  easily  described  explicitly.  First  of  all,  (7) 
shows  that  the  velocity  field  of  any  such  solution  must  be  horizontal, 
i.e.  directed  along  surfaces  of  constant  9,  u  •  VCp  =  0.  (We  shall 
always  use  the  word  "horizontal"  in  this  sense  -  vectors  perpendicular 
to  k  are  thus  not  usually  horizontal.)  If  this  is  true,  then  the 
cross  product  of  V<p  with  the  time  independent  form  of  (5)  shows 
that 

u  =  (2CpzP0)_1Vq>  X  Vp  (12) 

and  putting  this  back  in  (5)  we  find  that  P  must  be  given  by 


These  two  necessary  conditions  together  with  the  equation  of  state 

(8)  give  us  a  candidate  for  a  steady  solution  of  (5)  -  (8)  for  any 

function  p,  and  in  fact  one  easily  checks  that  formulas  (12) 

and  (l],)  are  sufficient  for  (5)  and  (7)  to  be  satisfied.  However 

(6)  will  be  satisfied  only  if  also  Vp  •  V9  xVCp  =  0;  this  is 

z 

no  restriction  on  p  in  case  Vcp  x  V9  =  0,  but  otherwise  it 

z 

shows  that  p  must  be  constant  along  the  integral  curves  of  the 

vector  field  V9  x  V9  ,  i.e.  along  the  intersections  of  the  surfaces 

z 

of  constant  9  and  those  of  constant  9  .  Parts  of  the  flow  region 

z 

R  throughout  which  V9  x  V9  =0,  i.e.  in  which  9  and  9  are 

z  z 

functionally  dependent,  we  shall  call  geostrophically  free  regions. 
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Elsewhere,  we  imagine  the  integral  curves  of  VcpxVcp  to  be  constructed, 

z 

and  call  them  geostrophic  curves.  Through  each  point  in  parts  of  R 

which  are  not  geostrophically  free  (if  any)  there  passes  a  unique 

geostrophic  curve;  some  or  all  of  these  curves  may  cross  the  boundary 

of  the  region  R.  Those  parts  of  R  which  are  covered  by  geostrophic 

curves  which  cross  the  boundary  B  we  call  geostrophically  blocked 

regions.  Other  parts  of  R  which  are  covered  by  geostrophic  curves 

which  do  not  cross  the  boundary  we  call  geostrophically  guided  regions. 

We  shell  assume  that  the  geostrophic  curves  in  guided  regions  are 

closed,  neglecting  the  possibility  that  they  might  wind  around  infinitely 

inside  the  region.  It  is  certainly  difficult  to  imagine  a  potential 

function  which  is  at  all  reasonable  for  which  the  geostrophic  curves 

stay  in  the  flow  domain  but  are  not  closed.  With  this  terminology 

we  may  say  then  that  for  a  steady  flow  it  is  necessary  that  p  be 

constant  along  geostrophic  curves  in  guided  or  blocked  regions,  i.e. 

(essentially)  p  =  p(<P,cp  )  there.  There  is  an  additional  minor 

z 

restriction  on  p  -  it  must  be  such  that  (12)  and  (13)  do  not  make 

u  or  p  singular  at  places  where  <P  =  0,  if  there  are  any  such. 

-  z 

But  except  for  this,  (12),  (13)  and  the  equation  of  state  (to  give 
s)  provide  the  general  solution  of  the  steady  forms  of  (5)  -  (8), 
with  p  arbitrary  in  free  regions,  or  an  arbitrary  function  constant 
on  geostrophic  curves  otherwise.  However  we  are  actually  interested 
only  in  solutions  which  also  satisfy  the  boundary  condition  u  •  n  =  0 
on  B;  (12)  shows  that  this  is  always  true  on  horizontal  parts  of 
B,  if  any,  but  that  it  is  necessary  (and  sufficient)  that  p  be 
constant  along  horizontal  curves  on  the  non -horizontal  part  of  B. 
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Such  a  steady  solution  of  (5)  -  (8)  which  satisfies  u  •  n  =  0  on  B 
we  call  a  geostrophic  flow.  In  non-free  regions  we  note  that 

!  &p(<P,<P  ) 

u.(2fl>2P0)  — 

z 

so  that  the  velocity  vector  of  a  geostrophic  flow  is  directed  along 

the  geostrophic  curves.  Furthermore  the  magnitude  of  u  dong  a 

given  geostrophic  curve  (in  the  non-free  case)  is  a  constant  multiple 

of  |vqp  xVqp  |,  and  thus  if  a  geostrophic  curve  crosses  the  boundary 
z 

the  velocity  vector  on  it  must  be  zero.  In  this  case  the  momentum 
equation  (or  (12))  shows  that  p  must  not  only  be  constant  along 
geostrophic  curves,  but  constant  on  horizontal  surfaces,  i.e. 
(essentially)  p  =  p(<p)  in  blocked  regions;  and  (15)  and  (8)  then 
show  that  p  and  s  are  likewise  constant  on  horizontal  surfaces. 
Thus  a  geostrophic  ‘flow"  in  a  blocked  region  is  at  most  only  a 
perturbation  of  the  basic  state  into  another  state  of  hydrostatic 
equilibrium. 

Summarizing  these  results,  the  geostrophic  flows  are  given  by 
(12),  (15)  and  (8)  with  p  a  function  which  is:  (a)  arbitrary  in 
free  regions,  except  that  it  must  be  constant  along  horizontal 
curves  on  the  non -horizontal  part  of  B,  (b)  an  arbitrary  function 
constant  along  geostrophic  curves  in  guided  regions,  or  (c)  an 
arbitrary  function  constant  on  horizontal  surfaces  in  blocked  regions. 
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4.  Conservation  theorems.  The  temporal  conservation  of  the 

energy  integral  is  a  result  of  a  very  general  sort;  something  like 

it  is  true  of  any  nou-dissipative  system.  We  come  now  to  some  results 

which  are  in  a  sense  much  stronger,  and  relate  much  more  specifically 

to  the  properties  of  linearized  rotating  stratified  flow.  Ihe  first 

of  these  has  to  do  with  quantities  on  the  boundary  B  of  the  flow 

region.  Some  of  B  may  consist  of  pieces  of  constant  <P  surfaces  - 

we  call  this  the  horizontal  part  B^  of  B;  on  it  n  x  Vcp  =  o. 

Parts  of  the  boundary  on  which  n  xVf/o  make  up  the  non-horizontal 

part  B^.  (Edges,  where  Bn  and  B^  meet,  we  may  count  as  part 

of  B^.)  Since  n  and  VCp  have  the  same  direction  on  the  "flat 

spots"  Bh,  the  boundary  condition  and  the  entropy  equation  (7) 

show  that  =  0  there.  Thus  entropy  is  conserved  pointwise  on 

flat  spots,  for  any  solution  satisfying  the  boundary  condition. 

This  is  not  true  on  B  ,  but  there  is  a  related  result.  Let  T 

n  n 

be  any  closed  horizontal  curve  on  B  ,  and  consider  the  rate  of 

n 

change  of  the  circulation  around  f  .  We  have 


St  Jr  a  •  «£- 

n 


[2k  x  u  +  p^Vp  +  PQlpVcp^  '  d- 


n 


=  -2  /  k  x  u  •  dr 

r 

n 


since  is  horizontal  (so  vcp  .  dr  =  0),  is  constant  along 

T  ,  and  Vp  •  dr  is  an  exact  differential.  However,  on  r  , 
nJ  *  -  n 


vcp  X  n 


dr  = 


-  Vcp  x  n 


da, 


da  being  arc  length,  so 
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n  •  k  n  •  k  s 

k  x  u  •  dr  =  -  -  u  •  VCp  da  =  - — -  da  • 

|vcp  x  n|  |v<p  x  n|  S0 


Thus  we  obtain 


n  •  k 
|vcp  x  n| 


so  setting 


C(T  )  =  /  [u  x  Vip  +  2  4  kl  •  H  — — —  (1^) 

n  so  Ivtpxnl 

n  - 

we  see  that  C(F  )  is  constant  in  time.  We  have  established  our 
'  n 

first  basic  result: 

Theorem  I.  For  every  solution  of  (5)  -  (8)  with  u  •  n  =  0 
on  B, 

(a)  s  is  conserved  on  the  flat  spots  B^,  if  any,  and 

(b)  C(r  )  (given  by  (I1*))  is  conserved  for  each  closed  horizontal 

curve  T  lying  on  B  . 
n  n 

Our  second  basic  conservation  property  concerns  a  quantity  related 
to  that  appearing  in  the  integrand  of  (14).  One  readily  shows  from 

(5)  and  (7)  that 

It  {%*?  -  X  VCP  +  ^O^  }  =  "2p0Ji  '  X  VCP*  (15) 
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Taking  the  divergence  of  this  and  using  (6)  we  get 


x  V<p  +  2s(s^) 


*  if v* 


Vcp  x  vcp  .  (16) 

z 


Let  us  define  the  scalar  field  n  by 


n  =  p-xv  •  [p^1  (uxV<p  +  2(s£) _1sk) 


-  2polp- 


(17) 


We  call  II  the  "potential  vorticity"  of  the  flow.  Throughout  any 
geostrophically  free  region  we  see  at  once  from  (l6)  that  II  is 
conserved  pointwise.  In  non -free  regions  this  is  not  true,  but  a 
somewhat  weaker  kind  of  conservation  does  apply.  First  consider 
a  geostrophically  guided  region,  and  let  r  be  any  closed  geostrophic 
curve  in  it  (or  on  its  boundary).  We  define  the  "mean  potential 
vorticity  for  r"  by  the  formula: 


n(r)  =  /  cp2 1  vcp  x  vcp  |_;Ln  do. 
jr  z  z 


(18) 


Along  T  we  have  dr  =  (vcp  x  vcp  |  Vcp  x  vcp  do,  so  (l6)  shows  that 


/  p”^  |  vcp  x  vcp  |  ”Vp  •  Vcp  x  vcp  do 
\J  -p  u  z  z 


p  1Vp  •  dr  =  0, 


since  pQ  is  constant  along  r  and  Vp  •  dr  is  exact.  This  kind 
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of  a  result  does  not  apply  in  blocked  regions,  since  there  the 
geostrophic  curves  are  not  closed.  For  a  blocked  region  we  consider 
any  horizontal  surface  E  lying  in  the  blocked  region,  and  let  T 
be  the  horizontal  curve  on  B  which  is  the  boundary  of  E;  E  and 
r  are  characterized  by  some  value  of  the  potential,  say  cp^.  Then 
we  construct  the  quantify 


where  the  orientations  of  r  and  E  are  related  as  usual  in  Stokes’ 
theorem,  i.e.  if  |v<p|_1vcp  is  regarded  as  the  unit  normal  to  E, 
on  T  we  have  dr  =  |  VCP  x  n|"*'Vcp  x  n  do.  Now  we  have  on  T; 


(u^  x  V<p  +  2st(s^)”1k)  •  n|vcp  x  n|-1( 


do 


(-2k  Xu-  p”Vp)  x  Vcp  -  2u  •  Vcpk 


n|  Vcp  x  n  ^do 


2cp  u  -  p’Vp  x  V<p  .  n|  Vcp  x  n|  -1do  =  -p'Vp  •  dr,  so 

Z  ”  u  "*  ***  *  *  *** 


•  VC?  *  V<^  I  vcp|  _1dS  + 

Q.X-  yj  ^  0  Z  Z  p 


f 

p  z  0 


p  •  dr 


-P^V  x  (cp’Vp)  •  vcp  I  vcp|_1ds  +  f  cp_1p  Vp  •  dr  =  0 

Z  v 


by  Stokes'  Theorem,  since  p^  is  constant  on  E.  This  finishes 


the  proof  of 
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Theorem  II.  For  every  solution  of  (5)  -  (8)  with  u  •  n  =  0 
on  B, 

(a)  II  (given  by  (17))  is  conserved  at  each  point  of  a  geostrophically 
free  region, 

(b)  II  (given  by  (18))  is  conserved  for  each  geostrophic  curve 
in  a  guided  region,  and 

(c)  Q  (given  by  (19))  is  conserved  for  each  horizontal  surface 
in  a  blocked  region. 

5.*  The  initial -value  problem.  Suppose  u(t),  p(t),  s(t), 

p(t)  is  the  solution  to  the  initial  value  problem  for  some  given 

set  of  initial  conditions.  We  shall  define  the  geostrophic  part 

of  u,  etc.,  (called  u  ,  etc.)  to  be  the  projection,  orthogonal 

6 

in  the  sense  of  the  inner  product  (11),  of  u,  etc.  onto  the  space 

of  geostrophic  flows.  Of  course  we  really  have  such  a  projection 

for  each  value  of  t,  but  these  are  in  fact  always  the  same . 

This  would  be  immediately  obvious  if  we  knew  the  existence  and 

completeness  of  normal  mode  solutions  in  terms  of  which  we  could 

represent  u(t)  etc.,  for  all  oscillatory  normal  modes  are 

orthogonal  to  the  steady  ones,  which  are  just  the  geostrophic 

flows.  But  even  without  the  exlscence  of  a  complete  set  of  normal 

modes,  the  projection  is  easily  seen  to  be  independent  of  time; 

this  follows  fo1'  instance  from  the  equivalent  definition  of  u^ 

etc.  as  that  geostrophic  flow  which  minimizes  the  energy  of  the 

difference  flow  u  -  u  ,  etc.,  (at  a  given  time  t).  If  u 

“8  6^ 

is  the  projection  at  t1  and  u  that  at  t  then  the  energies 
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of  the  difference  flows  u  -  u  and  u  -  u  must  be  constant 

“  "Bl  -  “g2 

in  time,  hence  equal ,  since  one  is  minimal  at  t^  and  the  other 

at  t».  The  identity  of  u  and  u  then  follows  from  the 
2  “g2 

uniqueness  of  the  projection  implied  by  the  positive  definiteness 
of  the  energy.  Since  the  geostrophic  part  is  independent  of  time, 
it  can,  in  principle  at  least,  be  determined  from  the  initial  data  - 
it  is  our  main  purpose  to  describe  explicitly  how  this  may  be  done. 

This  problem  is  not  altogether  trivial  because  the  space  of  geostrophic 
flows  is  relatively  extensive  -  it  is  not  simply  a  matter  of  calculating 
a  coefficient,  as  it  would  be  if  there  were  only  one  geostrophic 
flow. 

It  is  a  reasonable  conjecture  that  the  geostrophic  part  might 

also  be  the  time  average  over  0  <  t  <  °°  of  the  solution  to  the 

initial  value  problem.  This  is  probably  essentially  true;  for 

example  if  one  assumes  that  u(t)  etc.  remain  bounded  and  that 

the  time  average  exists,  one  easily  checks  that  the  time  average 

satisfies  the  geostrophic  equations.  Then  the  time  average  of 

u  -  u  etc.  would  also  be  a  geostrophic  flow,  and  one  that  is 
—  6 

orthogonal  to  all  geostrophic  flows  (assuming  interchangability  of 
spacial  and  temporal  integrations),  in  particular  to  itself.  Thus 
the  time  average  of  u  would  in  fact  be  u^.  It  should  be  remarked 
however  that  this  argument,  even  if  it  could  be  fully  justified, 
would  only  prove  that  u  was  the  time  average  almost  everywhere. 

It  seems  in  fact  to  be  the  case  that  sometimes  the  time  average  is 

not  exactly  the  same  as  u  ,  but  exhibits  certain  discontinuities 

"g 

on  parts  of  the  boundary.  We  shall  see  later  how  this  comes  about. 
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To  determine  u  we  first  need  a  formula  for  the  inner  product 
“g 

of  an  arbitrary  time  dependent  flow  u,  p,  p,  s  with  an  arbitrary 
geostrophic  flow  u,  p,  p,  s;  the  latter  is  given  in  terms  of  its 
pressure  field  p  by  the  formulas  of  §3,  and  we  want  to  express 

A 

the  inner  product  also  in  terms  of  p.  We  have,  using  the  real  form 
of  (11)  (the  complex  form  is  only  needed  in  connection  with  oscillatory 
normal  modes) : 


(u,p,p,s|u,p,p,s)  =  J  PqU  •  u+  P^fpPP  -  P~1(s(^)‘1ssJ  dV 

-  JR  [l  ^  x  •  a  ♦  Polfp»  +  \\X  *  v)]  dv 

"  IR  k1  (  2  a  *  vtp  +  (si5 ’S)  ‘  Vp  +  V  (polp  +  (si>'ls). 

u  x  vcp  +  2(Sq) _1sk  dS  +  J  fpp  +  (s^)'V)  dV 
qT1  (u  x  V<p  +  2(s^)'1sk)  dV. 


1  m"l^ 

-  cpz  pn  . 


|pv 


But 


■p;1  (a  x  V<p  +  2(s1J)"1sk)j  -  Pg1? 


LPolpz1  (a  xVtp  +  2(sorh± 


•  2c0lpi<,,)(si)'ls  ■  "  +  2p0lp  •  2P01(fsSi  -  Vo,(sirls 


=  n  +  2polfp(p  +  po(so>’ls) 


Thus  we  obtain: 
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pTE  dV  +  /  9~1p(s')~1sn  •  MS 

J  tj  z  U 


\ 


q>_1pn  •  [u  x  Vcp  +  gCsM^skjdS. 

Z  *  0 


(20) 


Now  p  is  the  pressure  field  of  a  geostrophic  flow,  so  it  must  be 
constant  along  the  horizontal  curves  r^(cp)  on  B^.  We  imagine 
the  last  surface  integral  in  (20)  to  be  calculated  by  dividing 
Bn  up  into  infinitesimal  strips  between  adjacent  curves  1^;  one 
easily  sees  from  the  geometry  that  in  this  case  we  have 
dS  =  |  Vcp  x  n|  ^d<P  do,  so 


/  <p-1pn  •  [u  xvcp  +  2(s')-1sk]dS 

Jb  Z  —  U 


n 


=  f  p(cp)d<p  [  CP-1  i V7<p  x  n|_1n  •  (u  x  Vcp  +  2(s » ) _1sk)do  (21) 

°  T  (<p)  z  U 

n  n 


We  now  use  these  formulas  to  characterize  those  flows  which  are 
orthogonal  to  all  geostrophic  flows,  by  requiring  the  right  hand  side 

/N 

of  (20)  to  be  zero  for  an  arbitrary  geostrophic  pressure  field  p. 

Let  us  suppose  the  region  3  to  be  divided  into  geostrophically 
free,  guided,  and  blocked  parts,  R^,  and  Rg,  and  denote  the 
portions  of  B^  and  Bn  which  are  associated  with  each  such  part 

/s 

by  similar  subscripts.  First  we  may  choose  p  to  be  zero  on  the 
boundary  and  in  R_  and  R_,  but  otherwise  arbitrary  in 

u  13 

(20)  then  shows  that  for  u,  etc.,  to  be  orthogonal  to  all 
geostrophic  flows  it  must  have  zero  potential  vorticity  throughout 
Rp.  Next  we  may  choose  a  p  which  is  zero  in  R^  and  Rg  (as  we 
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have  just  seen  the  volume  integral  then  vanishes  since  II  must 
be  zero  in  R^)  and  which  is  zero  on  Bn,  and  B^,  but 

which  is  arbitrary  on  B^.  We  thus  see  that  it  is  also  necessary 

that  s  should  vanish  (pointwise)  on  B.  _.  Now  we  cannot  choose  a 

Hr 

geostrophic  pressure  field  which  is  arbitrary  on  B^,  for  there 

it  must  be  constant  along  horizontal  curves,  but  its  value  can  vary 

arbitrarily  from  one  horizontal  curve  to  the  next.  Using  the 

form  (21)  for  the  surface  integral  over  B  „  we  see  then,  noting 

nr 

that  <P  is  constant  along  horizontal  curves  on  B  that  it 
z  nF 

is  also  necessary  that  the  quantity  C  defined  by  (14)  should 

be  zero  on  every  horizontal  curve  on  B^.  Next  we  may  choose  a 

p  which  is  zero  on  the  boundary  and  ir  R  but  is  arbitrary  in 

15 

R„,  to  the  extent  that  is  possible,  i.e.  an  arbitrary  function 
of  geostrophic  curves  in  Rg.  Since  p  is  not  completely  arbitrary 
we  cannot  conclude  that  H  must  vanish  in  R  ,  but  if  we  think  of 
the  volume  integral  over  R„  as  being  evaluated  by  integrating  over 

VJ 

infinitesimal  tubes  around  the  geostrophic  curves,  such  a  tube  being 

given  for  instance  by  cp  <  cp  <  cp.  +  dcp  and  cp  <  cp  <  cp  ,  +  d<p  , 

i  —  -  1  Z.L  —  z  —  zl  z 

one  sees  that  dV  =  |vtp  x  Vcp  |  ^dcp  dcp  da,  and  the  vanishing  of 

z  z 

A 

the  volume  integral  for  arbitrary  5(9*9  )  thus  implies  that  the 
mean  potential  vorticity  TE( T)  defined  by  (18)  must  vanish  for 
every  closed  geostrophic  curve  in  R_.  As  above  one  then  shows  that 

U 

C(T  )  must  also  vanish  for  T  ’s  on  B  since  cp  is  constant 
v  n7  n  nG'  z 

along  these  curves,  the  geostrophic  curves  in  Rn  not  crossing  the 
boundary.  However  one  cannot  conclude  that  s  =  0  pointwise  on 

A 

B,  _,  for  p  must  be  constant  along  the  geostrophic  curves  there. 

nu 


.-fir  *■*? 
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Dividing  B^  into  strips  between  adjacent  geostrophic  curves  ve  find 
dS  =  |vtp  x  Vcp  |"^|vcpjd<P  da  and  thus  since  n  •  k  =  +  9  |vcp|  ^  on 
B^,  we  see  that  for  each  geostrophic  curve  on  the  quantity 


D(rh)=/  Iwx 


V<p 


i-l 


s  da 


(21) 


must  vanish.  Passing  now  to  the  blocked  parts,  if  any,  we  note  that 

A 

here  the  arbitrariness  of  p  is  so  restricted  that  we  can  no  longer 
consider  the  volume  integral  and  the  surface  integral  over 
separately.  Evaluating  the  former  by  integrating  over  horizontal 
sheets  between  adjacent  surfaces  of  constant  9  and  the  latter  by 
using  the  corresponding  horizontal,  strips  on  B  we  find  that  the 
necessary  condition  is  the  vanishing  of  the  quantity  Q(E)  defined 
by  (19),  for  each  horizontal  surface  £  in  R._.  Finally  on  B,^, 

d  n_D 

since  n  •  k  =  +  |v<p|  ^9  there,  we  must  have  the  quantity 
—  —  —  z 


E  = 


|  V9|"1s  dS 


(22) 


zero  as  well.  It  should  perhaps  be  mentioned  explicitly  that  in 
case  some  of  the  surfaces  £  or  B^  consist  of  several  disconnected 
pieces,  these  conditions  must  hold  separately  for  each  connected 
component.  Summarizing  these  results  we  have 


Theorem  III.  In  order  that  a  flow  u,  p,  P,  s  should  be 
orthogonal  to  all  geostrophic  flows  it  is  necessary  (and  sufficient) 


that  each  of  the  following  should  vanish: 
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(a)  The  potential  vorticity  II  (eq.  (17) ),  at  each  point 
of  Rj.; 

(b)  n(T)  (eq.  (18)),  for  each  geostrophic  curve  r  in  R„; 

Ci 

(c)  Q(5!)  (eq.  (19)),  for  each  horizontal  surface  in  R^; 

(d)  C(T  )  (eq.  (14)),  for  each  horizontal  curve  on  B  _ 

n  nr 

or  BnG* 

(e)  s,  on  B^; 

(f)  D(r^)  (eq.  (21)),  for  each  geostrophic  curve  on  B^q, 

(g)  E  (eq.  (22)),  for  each  connected  component  of  B  . 

xlD 

Notice  that  all  of  these  quantities  are  constant  in  time, 
because  of  Theorems  I  and  II.  They  must  of  course  be,  because  of 
the  temporal  invariance  of  the  inner  product  which  assures  us  that 
if  u  etc.  is  once  orthogonal  to  all  geostrophic  flows  it  will 
always  be  so;  note  however  that  Theorem  I  is  slightly  stronger 
since  it  asserts  the  constancy  of  s  pointwise  on  all  of  B^, 
without  the  averaging  implied  in  the  quantities  D  and  E. 

Theorem  III  now  allows  us  to  characterize  completely  the 
geostrophic  part  of  the  solution  to  the  initial  value  problem. 

For  u  etc.  to  be  the  projection  of  u  etc.  onto  the  geostrophic 

"8 

flows  simply  means  that  the  difference  flow  u  -  u  etc.  is  orthogonal 

"8 

to  all  geostrophic  flows,  hence  the  quantities  (a)  •••  (g)  mentioned 
in  Theorem  III  must  all  vanish  when  calculated  for  this  difference 
flow.  Thus  Ug  etc.  must  be  a  geostrophic  flow  for  which  the 
quantities  (a)  •••  (g)  are  all  the  same  as  they  are  for  the  initial 
data.  The  geostrophic  part  is  uniquely  determined  by  these  conditions, 
for  the  difference  of  two  geostrophic  flows  both  of  which  have  the 


■ .  ,  ov 


same  values  of  (a)  •••  (g)  would,  be  a  geostrophic  flow  with  zero 
values  of  these  quantities,  hence  orthogonal  to  itself  and  so  zero. 
This  completes 

Iheorem  IV.  The  geostrophic  part  of  the  solution  to  the  initial 
v-ni ne  problem  is  that  (unique)  geostrophic  flow  which  has  the  same 
values  of  Lhe  quantities  mentioned  in  Theorem  III  as  the  initial 
flow  has. 

The  actual  computation  of  u  may  not  be  easy,  but  it  is  easier 

6 

than  the  full  initial  value  problem.  For  example  if  the  entire  region 
is  geostrophically  free,  the  condition  II  =  H  gives  an  elliptic 
partial  differential  equation  in  3  variables  for  p  ,  for  which 
in  addition  we  have  the  boundary  conditions  that  pg  is  constant 
on  horizontal  curves  on  Bn,  C  =  CQ  on  these  curves,  and  s  =  sQ 
on  B^.  We  do  not  have  a  proof  of  the  existence  of  a  solution  to 
this  problem  (as  usual,  uniqueness  is  much  easier)  but  we  shall  see 
an  example  in  the  next  section  which  at  least  suggests  that  such 
existence  is  to  be  expected. 

The  unique  determination  of  the  geostrophic  part  by  the  initial 
values  of  the  quantities  (a)  •••  (g)  indicates  that  we  may  describe 
them  as  "the  features  of  the  initial  data  which  are  carried  by  the 
geostrophic  part."  It  is  natural  to  expect  that  the  geostrophic  part 
would  carry  all  aspects  of  the  initial  data  which  are  constant  in 
time  -  yet  this  is  not  quite  so,  at  least  when  guided  or  blocked 
regions  are  present,  for  Theorem  I  shows  for  instance  that  s  is 


1*0*  'r'y 
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pointwise  constant  on  B^,  yet  the  geostrophic  part  can  only  carry 
certain  averages  of  s  on  and  B^;  similarly  C(r^)  is 

conserved  for  *11  s  on  Bn,  but  the  geostrophic  part  may  not 
be  able  to  carry  all  of  C  on  B  though  it  does  on  the  rest 

HU 

of  Bn>  It  seems  however  that  this  somewhat  peculiar  situation 
occurs  only  at  the  boundary;  at  any  rate  simple  examples  which  can 
be  explicitly  solved  and  in  which  similar  phenomena  occur  indicate 
this,  and  we  have  seen  above  that  with  fairly  plausible  assumptions 
one  can  show  the  identity  of  the  geostrophic  part  and  the  time 
average  almost  everywhere. 


5..  Examples  and  remarks.  Let  us  consider  first  the  situation 
typical  of  most  laboratory  experiments  where  the  gravitational  field 
is  that  of  the  earth;  to  have  it  constant  in  the  rotating  system  we 
must  assume  that  the  vectors  H  and  £  are  aligned,  so  that  the 
dimensionless  potential  function  has  the  form 


2 

where  y  =  g/(H 1  ).  In  this  case  the  entire  region  is  geostrophically 

free.  (If  fi  and  g  are  not  parallel,  we  will  not  of  course  have 

any  steady  basic  state;  this  is  indeed  one  of  the  difficulties  in 

doing  experiments  with  rotating  stratified  fluids  -  the  alignment 

must  be  quite  precise  or  it  will  be  impossible  to  produce  a  sufficiently 

stationary  basic  state.)  Most  laboratory  experiments  are  also  done 

with  fluids  for  which  some  sort  of  Boussinesq  approximation  is  appropriate. 
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In  our  non-dissipative  case  we  may  take  this  to  mean  that  variations 

in  density  due  to  pressure  alone  are  negligible  (f  =  0;  this 

eliminates  compression  waves)  and  variations  in  density  due  to 

heating  (or  salt;  s  may  be  taken  to  be  salinity  rather  than 

entropy  if  the  stratification  is  produced  that  way)  are  small 

compared  to  the  average  density.  A  convenient  simple  model  is 

obtained  by  taking  the  basic  stratification  to  be  such  that  s^ 

is  constant,  and  assuming  that  also  f  =  -a  is  constant.  With 

s 

a  suitable  choice  of  the  basic  scales  of  density  and  entropy  we 
may  take  pQ  -  1,  s^  =  7  \  and  P  =  -Os  as  the  linearized 
equation  of  state,  p  may  then  be  eliminated  from  the  problem, 
and  the  basic  equations  of  this  model  become 

+  2k  x  U  +  Vp  -  ays(k  -  y~X(xi  +  y^))  =  0  (24) 

V  •  u  =  0  (25) 

st  +  u  •  (k  -  Y-1(xi  +  yj.))  =  0.  (26) 

This  may  be  still  further  simplified  by  taking  y  to  be  large, 
supposing  that  ay  =  N2  is  of  order  unity.  This  makes  the  horizontal 
surfaces  perpendicular  to  k.  In  this  case  the  energy  is 


(27) 
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and  the  potential  vorticity  is 


II  =  V  .  [u  x  k  +  2sk]  =  k  •  V  x  u  +  2s 


The  geostrophic  flows  are  given  by 


u  =  -  k  x  Vp 


„-2 

s  =  N  p„ 


for  any  function  p  constant  along  horizontal  curves  on  B^.  The 
potential  vorticity  of  a  geostrophic  flow  is  thus  given  by 


n  =  \  +  4n"2p  \ 

g  2  \  H  ^zz/ 


where 


V  =  V  -  k 
H  -  dz- 


The  quantity  C  of  (14)  is 


C(r  )  =  /  [u  x  k  +  2sk]  •  njk  x  n|  <3o 
n  j  p  ~ 


which  for  a  geostrophic  flow  becomes 


C  =  —  f  [V  p  +  4N~2p  k]  •  n|k  x  n|-1d0. 


2Jr  L  H' 
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A  specific  problem  for  which  explicit  calculations  can  readily  be 

made  is  obtained  with  this  model  by  talking  the  region  R  to  be  the 

2  2  2 

cylinder :  0  <  x  +y  <  a  ,  0  <  z  <  1.  To  compute  the  geostrophic 

pant  of  the  solution  to  the  initial  value  problem  we  compute 

IL  =  k  ■  V  x  u.  +  2s_  from  the  initial  data,  and  also  (going  to 
0  -0  Oz 

cylindrical  coordinates) 


r2ir 

CQ(z)  =  J  uQ(a,elZ)  •  ^  ad9. 


Hie  problem  for  the  geostrophic  pressure  field  p  is  thus 


+  Un'2i 


7-p  +  4N  p  =  211 
*zz  0 


with  boundary  conditions 


p  =  N^s-  on  z  =  0,1 
z  0 


p  =  const,  on  r  =  a  for  each  fixed  z 


(a,z,9)a  d9  =  2CQ(z) 


The  somewhat  unusual  boundary  conditions  (56)  and  (37)  on  the 
non-horizontal  boundary  become  familiar  ones  when  p  is  represented 
as  a  Fourier  series  in  9: 


\  (m)/  \  im9 

P  =  /  P  ;(r,z)e 
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(36)  thus  is 


P^(a,z)  =0  u  /  0 


and  (37)  gives 


d£ 


(0) 


3r  (a,z)  =  77-  C0(Z)‘ 


To  verify  the  existence  of  a  solution  to  the  Neumann  problem  for 
p^  it  is  only  necessary  to  check  that  the  source  strength  2II^0^ 
in  (3k)  is  consistent  with  the  net  flux  out  through  the  boundaries 
given  by  (35)  and  (37);  this  is  readily  done.  The  solution  then 
can  be  computed  for  example  by  expanding  p^  -  ra/m  CQ(z)  in 
a  series  of  the  form  £  a  (z)J„(k  r/a),  where  the  k  are  roots 
of  Jq,  and  p^  in  series  £  aj^  (z)  J|m|  (kj^r/a) ,  where  the 
k ^  are  roots  of  J|m|*  An  alternative  representation  is  obtainable 
by  removing  the  inhomogeneity  in  the  boundary  condition  (35)  instead 
of  that  in  (37),  and  using  Fourier  cosine  series  in  z. 

With  this  simple  model  it  is  in  fact  possible  to  find  explicitly 
all  the  oscillatory  normal  modes  as  well,  and  so  solve  the  full  initial 
value  problem.  Although  this  brings  out  a  number  of  interesting 
features,  we  do  not  do  this  here,  our  main  interest  in  this  paper 
being  in  the  geostrophic  part. 

Interesting  examples  which  can  be  solved  explicitly  to  the  same 
ex.ent  as  our  first  one  do  not  seem  to  be  so  easily  found  in  the 
geostrophically  guided  case,  but  without  attempting  a  full  discussion 
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we  may  consider  the  following.  We  take  the  same  sort  of  Boussinesq 
fluid  as  above,  but  now  the  region  will  be  a  spherical  shell  given 
in  dimensionless  form  (spherical  coordinates)  by  a  <  r  <  a  +  1. 

For  the  geopotential  we  shall  take  <P  =  -'ya  /r  so  that  "gravity" 
points  radially  inward,  and  as  above  we  assume  uniform  stratification, 
with  scaling  such  that  s^  =  y"1.  We  also  suppose  y  is  large  and 

a  small  so  that  N2  =  ay  is  of  order  unity,  lhen  we  find 

2  2 

Vcp  =  r^ya  /r  ,  where  r^  is  the  radial  unit  vector;  and 

<p  -  ya2z/r^.  Thus  V<p  x  Vcp  =  r,  x  ky^ aVr'*, 
z  '  z  —1  — 

and  the  geostrophic  curves  are  latitude  circles  on  the  spheres  of 
constant  r.  The  potential  vorticity  is  found  to  be 


0  being  the  polar  angle.  The  geostrophic  flows  are  given  by 

=  |  sec  6  ^  x  Vp  (39) 

and 

r~)  ^  .  j  ^  ^  o 

s^  =  r^a  “N  “sec  9  =  r“a  CN  ^(Pr  -  tan  0  Pg/r)  (40) 

where  p  =  p(r,0)  is  constant  on  the  geostrophic  curves,  but  is 
otherwise  arbitrary  since  the  boundaries  r  ---  a,  a  +  1  are  entirely 
horizontal.  Computing  the  potential  vorticity  for  a  geostrophic 


flow  one  obtains ; 
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Since  the  weight  factor  in  (18)  is  independent  of  azimuth, 

as  is  II^,  the  basic  equation  =  nQ  is  equivalent  to  setting 

IT  equal  to  the  average  over  azimuth  of  the  initial  potential 
8 

vorticity  field.  Similarly  the  boundary  conditions  on  p, 

2  -2  “2 

(cf.  (21)),  are  equivalent  to  taking  r  a  N  (p  -  tan  0  pQ/r)  to 
be  the  azimuthal  average  of  the  initial  entropy  on  r  =  a  and  a  +  1. 
The  explicit  calculation  of  the  geostrophic  part  of  the  flow  produced 
by  given  initial  conditions  appears  to  be  a  fairly  difficult  problem 
in  general  with  this  model,  but  by  carefully  choosing  the  initial 
conditions  it  is  reasonably  easy  to  give  some  simple  examples.  One 
class  of  such  examples  is  obtained  with  initial  conditions  of  the 
form: 

Uq  -  0,  s^  =  A(r)cos^S  +  B(r)cos  0 


(42) 
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(with  some  restrictions  on  A  and  B  to  be  mentioned  in  a  moment), 
or  any  initial  conditions  whose  azimuthal  averages  are  of  this  form. 
In  this  case  the  potential  vorticity  is  found  from  (38)  to  be 


nQ  =  2r^a  2A’(r)cos^9  +  2a  ^[(r^B)’  +  2rA]cos  0. 


A  little  calculation  shows  that  a  geostrophic  flow  determined  by 
p  =  f(r)cos^9  will  give  this  same  potential  vorticity  field  provided 
that  A  and  B  are  related  to  f  by 


A  =  a_2N'2r5(f/r5) 


B  =  3a”2N'2rf  -  (3a2/2r2)  f  t^r^r'2^ 


Such  a  geostrophic  flow  has  entropy 


s  =  a  2N  r^(f/r  )  'cos^O  +  3rf  cos  ©1 
6  L  J 


and  thus  as  far  as  the  cos^©  term  is  concerned  the  conditions  that 
s  =  s  on  r  -  a  and  a  +  1  are  satisfied  for  any  f ;  but  the 

6  vj 

cos  9  term  shows  that  we  also  must  require  that 


f(r)r  2dr  =  0. 
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Thus  to  construct  an  example  we  need  only  take  some  f  satisfying  this 
condition,  and  calculate  A  and  B  from  (45)  and  (46).  f  cos^9 
will  then  be  the  pressure  field  of  the  geostrophic  part  of  the  flow 
that  evolves  from  the  initial  state  (42).  Its  velocity  field  is, 
from  (39) » 


u  =  -  \  sin  29(f/r)2  (47) 

s  ^ 

where  <jP^  ^he  unl't  azimuthal  vector;  the  entropy  (46)  may  also 
be  written  in  the  form: 

sg  =  sQ  +  |  a2r“2cos  0  j  f(r1)r"2dr1.  (48) 

A  particularly  simple  explicit  example  is  obtained  by  taking 
f  =  r  (r  -  a  -  —  )/a  ,  in  this  case,  if  W  is  moderate  and  a 
is  rather  large  (thin  shell)  the  initial  entropy  field  is  not  too 
different  from  cos^©/(aN2) ,  and  is  thus  a  warming  of  the  northern 
and  cooling  of  the  southern  hemispheres.  The  geostrophic  velocity 
field  (47)  is  then  vertically  sheared,  being  from  west  to  east  at 
low  altitudes  and  from  east  to  west  at  high  altitudes  in  the  northern 
hemisphere,  and  reversed  in  the  southern. 

It  is  of  interest  to  compare  our  general  results  with  their 
analogues  for  a  rotating  homogeneous  fluid,  for  the  relationship 
between  the  two  is  a  little  more  subtle  than  one  might  perhaps 
naively  expect.  Of  course  one  would  anticipate  that  "turning  on" 
a  little  stratification  (with  a  gravitational  field)  would  be  a  small 
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perturbation,  at  least  for  a  finite  time,  on  the  full  initial  value 
problem,  but  this  is  not  necessarily  the  case  when  attention  is 
focussed  on  the  steady  part  of  the  solution..  Indeed  steady  flows  in 
the  linear  theory  of  rotating  homogeneous  fluids  must  be  independent 
of  z  ( Taylor -Proudman  Theorem),  and  this  is  not  .t  all  necessary 
for  steady  stratified  flows ;  on  the  other  hand  steady  stratified 
flows  must  be  along  geopotential  surfaces,  which  is  not  necessary 
in  the  homogeneous  case  where  the  geopotential  (in  the  absence  of 
free  surfaces)  plays  no  role  at  all.  What  is  happening  here  is 
of  course  that  some  of  the  oscillatory  normal  modes  in  the  stratified 
case  have  frequencies  which  tend  to  zero  as  the  stratification 
disappears  and  so  show  up  as  geostrophic  motions  when  the  stratification 
is  turned  off.  Likewise,  some  initial  conditions  which  would  excite 
time  dependent  motion  (probably  largely  at  low  frequencies)  in  the 
homogeneous  case,  may  in  fact  lead  to  steady  flows  when  a  little 
stratification  is  put  in.  Though  the  situation  is  really  more  complex, 
it  is  a  little  like  the  apparent  discontinuity  in  the  nature  of 
geostrophic  homogeneous  flows  when  the  bottom  of  a  container  of 
constant  height  is  slightly  deformed:  the  class  of  geostrophic  flows 
is  greatly  curtailed,  and  new  low  frequency  oscillatory  modes  (Rossby 
waves)  appear. 

The  homogeneous  analogues  of  our  results  can  of  course  be  obtained 
by  the  same  methods  -  the  difference  comes  only  in  the  first  step,  the 
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description  of  the  geostrophic  flows  themselves.  Although  these  results 
are  already  fairly  well  known  (they  are  essentially  given,  for  instance, 
in  H.  P.  Greenspan's  monograph  "The  theory  of  rotating  fluids," 

Cambridge,  1968)  it  is  perhaps  worth  recalling  them  briefly  to 
emphasize  that  the  analogy  with  the  stratified  case  is  in  fact 
quite  close.  To  describe  the  homogeneous  geostrophic  l'J ows  one  must 
also  distinguish  free,  guided,  and  blocked  regions.  Hie  free  regions 
are  those  of  constant  height  (measured  parallel  to  the  rotation  axis); 
the  guided  regions  are  those  not  of  constant  height,  but  in  which 
the  contours  of  constant  height  ("geostrophic  curves")  do  not  cross 
the  boundaries;  and  the  blocked  regions  are  those  in  which  they  do. 

The  geostrophic  flows  in  free  regions  are  independent  of  z  and 
parallel  to  the  top  and  bottom,  but  otherwise  arbitrary;  they  are 
essentially  characterized  by  the  vertical  vorticity  component,  a  quantity 
which,  when  vertically  averaged,  is  conserved  for  any  motion  in  a  free 
region.  The  geostrophic  flows  in  guided  regions  are  independent 
of  z  and  along  geostrophic  curves;  they  are  characterized  by  their 
"mean  circulation,"  the  vertically  averaged  circulation  around  the 
closed  geostrophic  contours,  another  conserved  quantity,  in  guided 
regions.  The  only  geostrophic  "flow"  in  a  blocked  region  is  zero. 

The  steady  part  of  the  solution  to  an  initial  value  problem  "carries" 
the  vertically  averaged  vertical  vorticity  of  the  initial  data  in  free 
regions  and  the  mean  circulation  of  the  initial  daca  in  guided  regions, 
just  as  in  the  stratified  case  it  carries  the  potential  vorticity 
and  the  mean  potential  vorticity.  In  over-all  structure  the  stratified 
geostrophic  flow  problem  is  really  very  similar  to  the  homogeneous  case, 
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even  though  the  latter  is  not  readily  derived  from  the  former  as  a 
limiting  case. 

What  is  the  relationship  between  our  function  II  said  the  familiar 
potential  vorticity  of  non-linear  shallow  water  theory?  The  latter, 
the  (absolute)  vorticity  divided  by  the  local  depth  of  the  fluid 
layer,  is  a  quantity  which  is  conserved  following  particles,  not 
conserved  at  a  particular  place  like  II  (or  II),  so  a  direct  analogy 
is  not  to  be  expected;  but  they  are  in  fact  closely  related.  The 
true  generalization  of  the  potential  vorticity  of  shallow  water  theory 
is  given  by  Ertel's  Theorem,  which  (in  a  non -rotating  system) 
asserts  the  constancy  following  particles  of  the  quantity  p  V  x  u  •  Vs. 
This  holds  for  an  inviscid  compressible  fluid  subject  to  a  conservative 
body  force  per  unit  mass,  with  s  any  scalar  field  (in  particular 
the  entropy,  if  there  is  no  heat  conduction)  which  is  itself  constant 
following  particles  and  in  addition,  if  Vp  x  Vp  ^  0,  is  a  function 
of  p  and  p.  Ertel's  Theorem  can  be  readily  proved  from  the 
Helmholtz  vorticity  theorem  (which  is  true  under  these  hypotheses) 
on  using  the  above  properties  of  s .  The  conservation  of  potential 
vorticity  in  shallow  water  theory  follows  from  Ertel's  Theorem,  since 
shallow  water  theory  is  mathematically  analogous  to  two  dimensional 
compressible  flow  of  an  isentropic  gas  with  7=2;  in  this  analogy 
the  'gas  density"  is  the  mass  per  unit  area  of  the  fluid  layer,  hence 
proportional,  to  the  local  depth.  The  gas  motion  being  two  dimensional 
and  isentropic,  we  may  take  s  =  z,  and  the  conservation  of  potential 
vorticity  follows.  Thus  the  quantity  p  V  x  u  •  Vs,  or  (in  a  rotating 
system)  p  ^(2Qk  +  V  x  u)  •  Vs  might  appropriately  be  called  "potential 
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vorticity"  for  compressible  stratified  flow.  In  our  linearized  theory 
this  is  given  in  dimensionless  form  by 


p’V  +  ep/P0)'1(2k  +  evx  u)  •  (s^vcp  +  evs) 


s 1  c  ^"cp 
0  0  yz 


2  +  e 


xu  ■  Vcp  +  2(s^)_1s2  -  2PQ1p^ 


m 


The  coefficient  of  e  in  this  expression  is  closely  related  to  H, 
and  our  basic  formula  ( l6 )  from  which  the  conservation  of  II  (or  n) 
was  derived  is  an  immediate  consequence  of  Ertel's  Theorem  in  the 
linearized  case.  Indeed  we  have  from  (17)  that 


n  =  cp  x  u  •  vcp  + 


2<so 


-l 


-  2pols 


+  v^:1) 

2 


u  x  Vcp  +  2p 


-1 

0 


s  I? 


(50) 


If  we  use  (4g)  and  (50)  and  linearize,  Ertel's  Theorem  says 


5~  n  -  vcp-1  .  n  X  vcp  -  20  s.  (p.cp^/s*) 
dt  z  -=t  0  t  dz  0  z  '  0 


+  2(sopoV^  •  v(sopo\}  =  °- 


Replacing  u,  and  s  in  this  formula  by  their  values  from  (5)  and 

■“u  X- 

(?)  one  recovers  almost  at  once  equation  (16) . 

The  function  n  is  thus  approximately  the  first  order  part  of 
the  "true"  potential  vorticity  p  V  x  u  •  Vs  -  it  has  been  defined 
somewhat  differently  in  order  to  obtain  a  quantity  which  is  constant 
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in  time  (in  the  geostrophicaUy  free  case)  at  a  fixed  place.  Although 
the  ,:true::  potential  vorticity  is  no  doubt  of  more  basic  physical 
significance,  being  in  fact  conserved  (following  particles)  even 


* 


in  the  non-linear  problem,  we  have  felt  it  appropriate  to  give  this 
name  to  the  function  II  which,  as  we  have  seen,  plays  a  fundamental 
role  in  the  linear  theory  of  rotating  stratified  flow. 
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